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Abstract 

We consider the problem of computing critical points of the restriction of a polynomial map to an 
algebraic variety. This is of first importance since the global minimum of such a map is reached at a 
critical point. Thus, these points appear naturally in non-convex polynomial optimization which occurs 
in a wide range of scientific applications (control theory, chemistry, economics,...). 

Critical points also play a central role in recent algorithms of effective real algebraic geometry. Ex- 
perimentally, it has been observed that Grobner basis algorithms are efficient to compute such points. 
Therefore, recent software based on the so-called Critical Point Method are built on Grobner bases en- 
gines. 

Let /i, . . . , /p be polynomials in Q[xi, . . . , Xn] of degree D, V C C" be their complex variety 
and TTi be the projection map {xi, . . . , x„) xi. The critical points of the restriction of tti to V are 
defined by the vanishing of /i , . . . , /p and some maximal minors of the Jacobian matrix associated to 
fi, . . . , fp. Such a system is algebraically structured: the ideal it generates is the sum of a determinantal 
ideal and the ideal generated by /i , . . . , /p. 

We provide the first complexity estimates on the computation of Grobner bases of such systems 
defining critical points. We prove that under genericity assumptions on /i, . . . , /p, the complexity is 
polynomial in the generic number of critical points, i.e. [D — . More particularly, in the 

quadratic case D — 2, the complexity of such a Grobner basis computation is polynomial in the number 
of variables n and exponential in p. We also give experimental evidence supporting these theoretical 
results. 



1 Introduction 

Motivations and problem statement. The local extrema of the restriction of a polynomial map to a real 
algebraic variety are reached at the critical points of the map under consideration. Hence, computing these 
critical points is of first importance for polynomial optimization which arises in a wide range of applications 
in engineering sciences (control theory, chemistry, economics, etc.). 

Computing critical points is also the cornerstone of algorithms for asymptotically optimal algorithms 
for polynomial system solving over the reals (singly exponential in the number of variables). Indeed, for 
computing sample points in each connected component of a semi-algebraic set, the algorithms based on 
the so-called critical point method rely on a reduction of the initial problem to polynomial optimization 
problems. In LlOlfTTl (see also |26l l27ll28l ), the best complexity bounds are obtained using infinitesimal 
deformation techniques of semi-algebraic geometry, nevertheless obtaining efficient implementations of 
these algorithms remains an issue. 

Tremendeous efforts have been made to obtain fast implementations relying on the critical point method 
(see ||35l [T5l [34l [32] |3T1 [161 [36l ). This is achieved with techniques based on algebraic elimination and 
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complex algebraic geometry. For instance, when the input polynomial system (F) : /i = • • • = /p = 
in Q[a:;i, . . . , a;„] satisfies genericity assumptions, one is led to compute the set of critical points of the 
restriction of the projection tti : (a;i, . . . , x„) xi to the algebraic variety V^(F) C C" defined by F; this 
set is denoted by crit(7ri, V^(F)). 

The set crit(7ri, V^(F)) is defined by F and the vanishing of the maximal minors of the truncated 
Jacobian matrix of F obtained by removing the partial derivatives with respect to xi. This system is highly- 
structured: algebraically, we are considering the sum of a determinantal ideal with the ideal (/i , . . . , /p) . 

In practice, we compute a rational parametrization of this set through Grobner bases computations 
which are fast in practice. We have observed that the behavior of Grobner bases on these systems does 
not coincide with the generic one. In the particular case of quadratic equations, it seems to be polynomial 
in n and exponential in p which meets the best complexity known bound for the quadratic minimization 
problem [9 25 1. Understanding the complexity of these computations is a first step towards the design of 
dedicated Grobner bases algorithms, so we focus on the following important open problems: 

(A) Can we provide complexity estimates for the computation of Grobner bases of ideals defined by such 
structured algebraic systems! 

(B) Is this computation polynomial in the generic number of critical points? 

(C) In the quadratic case, is this computation polynomial in the number of variables (and exponential in 
the codimension)? 

Under genericity assumptions, we actually provide affirmative answers to all these questions. 

Computational methodology and related complexity issues. Grobner bases are computed using 
multi-modular arithmetics and we will focus only on arithmetic complexity results; so we may consider 
systems defining critical points with coefficients not only in Q but also in a prime field. 

Let K be a field, K be its algebraic closure and F = (/i, . . . , fp) be a family of polynomials in 
K[a::i , . . . , Xn] of degree D and V{¥) be their set of common zeroes in K". 

We denote the Jacobian matrix 

dfi ... dfi 

dx\ dxn 

dfp . . . dfp 

dxi dxn 

by jac(F) and the submatrix obtained by removing the first i columns by jac(F, i). The set of maximal 
minors of a given rectangular matrix M will be denoted by MaxMinors(M). 

Finally, let I(F, 1) be the ideal (F) + (MaxMinors(jac(F, 1))). When F is a reduced regular sequence 
and V{F) is smooth, the algebraic variety associated to I(F, 1) is exactly crit(7ri, V(F)). 

So, to compute a rational parametrization of crit(7ri , V^(F)), we use the classical solving strategy which 
proceeds in two steps: 

(i) compute a Grobner basis for a grevlex ordering of I(F, 1) using the algorithm (see (TT)); 

(ii) use the FGLM algorithm |fT8l[T9l to obtain a Grobner basis of I(F, 1) for the lexicographical ordering 
or a rational parametrization of a/I(F, 1). 

Algorithm F5 (Step ( i)) computes Grobner bases by row-echelon form reductions of submatrices of the 
Macaulay matrix up to a given degree. This latter degree is called degree of regularity. When the input 
satisfies regularity properties, this complexity of this step can be analyzed by estimating the degree of 
regularity. 

FGLM algorithm |18| (Step (ii)) and its recent efficient variant 1191 are based on computations of 
characteristic polynomials of linear endomorphisms in ]K[a:;i, . . . , Xn]/ I(F, 1). This is done by performing 
linear algebra operations of size the degree o/I(F, 1) (which is the number of solutions counted with 
multiplicities). 

Thus, we are faced to the following problems: 

(1) estimate the degree of regularity of the ideal generated by the homogeneous components of highest 
degree of the set of generators F, MaxMinors(jac(F, 1)); 
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(2) show that the above estimation allows to bound the complexity of computing a grevlex Grobner basis 
ofI(F, 1); 

(3) provide sharp bounds on the degree of the ideal I(F, 1). 

As far as we know, no results are known for problems (1) and (2). Problem (3) has already been investigated 
in the literature: see |33 | where some bounds are given on the cardinality of crit(7ri, V^(F)). We give here 
a new algebraic proof of these bounds. 

Main results. Let K[xi, . . . , Xn]D denote {/ e K[xi, . . . , Xn] \ deg(/) = D} and remark that 
it is a finite-dimensional vector space. In the following, we solve the three aforementioned problems 
under a genericity assumption on F: we actually prove that there exists a non-empty Zariski open set 

C K[xi , . . . , x„]^ such that for all ¥ £ ff: 

(1) the degree of regularity of the ideal generated by the homogeneous components of largest degree of 
F, MaxMinors(jac(F, 1)) is d,.og = D{p -l) + {D-2)n + 2 (see Theorem[I); 

(2) with the ^5 algorithm, the highest degree reached during the computation is bounded by dreg (see 
TheoremlU; 

(3) the degree ofI(F,l) is <S ^ DP{D- 

The degree of regularity given in fi j is obtained thanks to an explicit formula for the Hilbert series of the 
homogeneous ideal under consideration (see Proposition [T}. This is obtained by taking into account the 
determinantal structure of some of the generators of the ideal we consider. The above estimates are the key 
results which enable us to provide positive answers to questions A, B and C under genericity assumptions. 

Before stating complexity results on the computation of critical points with Grobner bases, we need to 
introduce a standard notation. Let w be a real number such that a row echelon form of a n x n-matrix with 
entries in K is computed within 0(n") arithmetic operations in K. 

We prove that there exists a non-empty Zariski open set C K[xi, . . . , a;„]^ such that for all F e 
^nIK[xi,...,x„]P: 

(A) computing a grevlex Grobner basis of I(F, 1) can be done within O (^("^^'^°'^) ^ arithmetic opera- 
tions in K (see Theorem[3]l; 

(B) computing a rational parametrization of crit(7ri, V{F)) using Grobner bases can be done within 
O ((5'* °^") arithmetic operations in IK (see Corollary|5); 

(C) when D ^ 2 (quadratic case), a rational parametrization of crit(7ri , V{F)) using Grobner bases can 
be computed within O (^("Jp^)" + ''^'^'^^ arithmetic operations in K, this is polynomial in n 
and exponential in p (see Corollary |3]l. 

We also provide more accurate complexity results. The uniform complexity bound given for answer- 
ing question (B) is rather pessimistic. The exponent 4.03w being obtained after majorations which are 
not sharp; numerical experiments are given to support this (see Section |6]l. Moreover, under the above 
genericity assumption, we prove that, when p and D are fixed, computing a rational parametrization of 
crit(7ri, V{F)) using Grobner bases is done within 0{D^-^'^") arithmetic operations in K (see Corollary 

s. 

We also give timings for computing grevlex and lex Grobner bases of I(F, 1) with the MAGMA com- 
putational algebra system and with the FGb library when K = GF(65521). These experiments show that 
the theoretical bounds on the degree of regularity and on the degree of I(F, 1) (Theorem|2]i are sharp. They 
also provide some indication on the size of problems that can be tackled in practice: e.g. when D ^ 2 and 
p — 3 (resp. D = 3 and p — 1), random dense systems with n < 21 (resp. n < 14) can be tackled (see 
Section|6]l. 

Related works. As far as we know, dedicated complexity analysis of Grobner bases on ideals defining 
critical points has not been investigated before. However, as we already mentioned, the determinantal 
structure of the system defining crit(7ri , V{F)) plays a central role in this paper 
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In II20II . we provided complexity estimates for the computation of Grobner bases of ideals generated by 
minors of a linear matrix. This is generalized in ^21] for matrices with entries of degree D. Nevertheless, 
the analysis which is done here differs significantly from these previous works. Indeed, in [2(1 2U ^ gener- 
icity assumption is done on the entries of the considered matrix. We cannot follow the same reasonings 
since MaxMinors(jac(F, 1)) depends on F. Nevertheless, it is worthwhile to note that, as in ||201|2T1 . we 
use properties of determinantal ideals given in [12]. 

Bounds on the number of critical points (under genericity assumptions) are given in ll33l using the 
Giambelli-Thom-Porteous degree bounds on determinantal varieties (see flT, Ex. 14.4.14]). 

In [9j, the first polynomial time algorithms in n for deciding emptiness of a quadratic system of equa- 
tions over the reals is given. Further complexity results in the quadratic case for effective real algebraic 
geometry have been given in [25]. In the general case, algorithms based on the so-called critical point 
method are given in [10, 11, 26 , 27,, 28] . Critical points defined by systems F, MaxMinors(jac(F, 1)) are 
computed in algorithms given in (|3] E] |5] IH |6l [35] [T] [161. Th^ RAGlib maple package implements the 
algorithms given in |[35l [T6l using Grobner bases. 

The systems F, MaxMinors(jac(F, 1)) define polar varieties: indeed, this notion coincides with critical 
points in the regular case). In ||3]|2] [5] |4] |6], rational parametrizations are obtained using the geometric 
resolution algorithm f23] and a local description of these polar varieties. This leads to algorithms comput- 
ing critical points running in probabilistic time polynomial in D''{p{D — 1))"^*'. Note that this bound for 
D = 2 and p = n/2 is not satisfactory. In this paper, we also provide complexity estimations for comput- 
ing critical points but using Grobner bases, which is the engine we use in practice. Our results provide an 
explanation to the good practical behavior we have observed. 

Organization of tlie paper. Section [2] recalls well-known properties of generic polynomial systems. 
Problems (1) and (2) mentioned above are respectively tackled in Sections [3] and |4] Problem (3) is solved 
at the end of Section [4] Complexity results are derived in Section [5] Experimental results supporting the 
theoretical results are given in Section |6] 

Conclusions and Perspectives. We give new bounds on the degree of regularity and an explicit formula 
for the Hilbert series of the ideal vanishing on the critical points under genericity assumptions. This leads 
to new complexity bounds for computing Grobner bases of these ideals. 

However, we only considered the unmixed case: all polynomials fi, . . . , fp share the same degree D. 
The mixed case (when the degrees of the polynomials fi, . . . , fp are different) cannot be treated similarly 
since the difference of the degrees induce a combinatorial structure which has to be investigated. We 
intend to investigate this question in future works using the Eagon-Northcott complex, which yields a 
free resolution of the ideal generated by the maximal minors of a polynomial matrix under genericity 
assumptions. From this, we also expect to obtain a variant of the F5 algorithm dedicated to these ideals. 

Acknowledgments. This work was supported in part by the GeoLMI grant (ANR 20 11 BS03 Oil 06) 
and by the EXACTA grant (ANR-09-BLAN-0371-01) of the French National Research Agency. 

2 Preliminaries 

Notations 1. The set of variables {xi, . . . , Xn} is denoted by X. For d (E N, Monomials((i) denotes the 
set of monomials of degree d in the polynomial ring IK[X] (where K is afield, its algebraic closure being 

denotedbyK). We let a denote the finite set of parameters {am : I < i < p,m e lJo<d<D Monomials(d)}. 
We also introduce the following generic systems: 

• -S — ifiT ■ ■ ■ ,fp) 'K{a)[X]P is the generic polynomial system of degree D: 

m monomial 
deg(m)<_D 

• = (fi , ■ ■ ■ ,fp) G K(a)[X]P is the generic homogeneous polynomial system of degree D: 

m monomial 
deg(m)— -D 



4 



We let V{F) C IK denote the variety ofF — (/i, . . . , fp). The projective variety of a homogeneous 
family of polynomials F'* is denoted by W{F^). The projection on the first coordinate is denoted by tti, and 
the critical points of the restriction of ni to V{F) are denoted by crit(7ri, V^(F)) C V{F). Also, I(F, 1) 
denotes the ideal generated by F and by the maximal minors of the truncated Jacobian TOafr« jac(F, 1). 

Throughout the paper, if R is a ring and I <Z R is an ideal, we call dimension of I the Krull dimension 
of the quotient ring R/ 1. 

The goal of this section is to prove that the ideal 1) is 0-dimensional. This will be done in 
Lemma 12] below; to do that we will use geometric statements of Sard's theorem which require K to have 
characteristic 0. This latter assumption can be weakened using algebraic equivalents of Sard's Theorem 
(see L14. Corollary 16.23]). 

Lemma 1. Let 1(5, 0) be the ideal generated by 5 and by the maximal minors of its Jacobian matrix. Then 
its variety 0)) C K(o) is empty and hence 1^(3") is smooth. 



Proof. To simplify notations hereafter, we denote by 
by removing their respective constant terms a^^'' , . . . 
rameters in /ii, . . . , /ip. Let denote the mapping 



hi, . . . ,hp the polynomials obtained from fi, . . . , fp 
, a^^' . We will also denote by ^ the remaining pa- 



c I — > {hi{c),...,hp{c)) 



Suppose first that ?/;(K(j2/) ) is not dense (for the Zariski topology) in K(ja^) . Since the image 
4>(K{s^) ) is a constructible set, it is contained in a proper Zariski closed subset 'W C ]K(^) . Since 
there is no algebraic relation between a^^\ . . . , a'f'^ and the parameters in £/, this implies that the variety 
defined by hi + aj^"* ~ ■ ■ ■ = hp + a^^"* is empty and consequently smooth. Since hi + a.^^ = fi, our 
statement follows. 

Suppose now that ip{K{£/) ) is dense in K(^) . Let Kq C K(j2/) be the set of critical values of ip. 
By Sard Theorem ['37' Chap. 2, Sec. 6.2, Thm 2], Kq is contained in a proper closed subset of K(j2/) . 
Again, there is no algebraic relation between oJ^\ . . . , a^^"* and the parameters in .s^. Consequently, the 
variety associated to the ideal generated by the system fi, . . . , fp and by the maximal minors of jac(5^) is 
empty. □ 

Corollary 1. Let 0) be the ideal generated by and by the maximal minors of its Jacobian matrix. 
Then the associated projective variety W{1{^^, 0)) C P"^^]K(a) is empty. 

Proof For 1 < i < n, we denote by Oi the set 

{(ci : . . . : c„) I c, ^ 0} C P"-1K(^ 
and we consider the canonical open covering of P"^^K(o): 

l<i<n 

Therefore W{I{d^', 0)) = Ui<j<„(T^(I(S''\ 0)) n O^). Denote by the system obtained by substituting 
the variable Xi by 1 in "S^. According to Lemma[T]applied to ^i, the variety V{I{^i, 0)) is empty. Therefore, 
the set W{I{d^, 0)) n is also empty. Consequently, 0)) =0. □ 

We can now deduce the following result. 

Lemma 2. The projective variety W{I{^'^, 1)) C P"~^]K(o) is empty, and hence dim(I(5''^, 1)) = 0. 
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Proof. We let (^o and ipi denote the two following morphisms: 



K(a)[a;i, . . 


■ 1 •^n\ 


K(a)[a;2,. 




g{xi,.. 


• •) '^n) 


h^- g(0,a;2,.. 


• ; '^n ) 


K(a)[xi, . . 


■ 1 


K(a)[a;2,- 






• 7 '^n) 


iH- g(l,a;2,.. 


• 7 -^n ) 



Then 1)) can be identified with the disjoint union of the variety 1))) C K{a) and 

the projective variety Wiipoili^^, 1))) C ¥"-'^Kia). 

• Notice that = o). Therefore, the ideal 1)) C K(a)[x2, . . . , x„] is 
spanned by (^i(F'') (which is a generic system of degree D inn — I variables) and by the maximal 
minors of its Jacobian matrix. According to Lemma[T] the variety V{(pi 1))) is empty. 

• Similarly, (/7o(I(5''^, 1)) — I{(po{'S^^),Q) C IK(a)[a;2, a;„] is generated by the homogeneous poly- 
nomials (fioi'S^) and by the maximal minors of the Jacobian matrix ia,c{(po{^^)). Thus, according to 
CoroUary[Tl the variety W{ipo{'i{d^ , 1))) is also empty. 

□ 

3 The homogeneous case 

In this section, our goal is to estimate the degree of regularity of the ideal 1(5^'^, 1) C K(o)[X] which is 
a homogeneous ideal generated by 5^'' and MaxMinors(5^'^, 1) (see Notations [T]l. Recall that the degree 
of regularity drcg(/) of a 0-dimensional homogeneous ideal / is the smallest positive integer such that 
all monomials of degree drcg(^) are in /. Notice that drcg(/) is an upper bound on the degrees of the 
polynomials in a minimal Grobner basis of / with respect to the grevlex ordering. 

Theorem 1. The degree of regularity of the ideal 1(5^^, 1) is 

d„g(I(^?^ 1)) = D{p ~1) + {D- 2)n + 2. 
Notations 2. To prove Theorem\l\ we need to introduce a few more objects and notations. 

• A set of new variables {uij '■ 1 < i < p,2 < j < n} which is denoted by U; 

• the determinantal ideal V C K[U] generated by the maximal minors of the matrix 

"1,2 • • ■ Wl,n 

• Qi, . . . , flp(„_i) G K(o) [U, X] which denote the polynomials Uij for 1 < i < p,2 < j < n 

aniiflp(„_]^)_l_]^, . . . , which denote the polynomials f'l, . . . , f^; 

• the ideals = 2? + (g^, . . . , g^) C K(o)[?7, X]; 

' if g K.[X] (resp. I C IK[X]j is a polynomial and -< is a monomial ordering (see e.g. I[13\ Ch. 2, 
§2, Def. 1]), LM^{g) (resp. LM^(/)j denotes its leading monomial (resp. the ideal generated by 
the leading monomials of the polynomials in I); 

• a degree ordering is a monomial ordering -< such that for all pair of monomials mi, m2 € IK[X], 
deg(mi) < deg(TO2) implies mi -< m2. 
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Obviously the polynomials g^, for 1 < fc < p{n — 1) will be used to mimic the process of substituting 

the new variables Uij by indeed we have 3(p„) n]K[X] — 1(3^ , 1). 

Our strategy to prove Theorem [Tlwill be to deduce the degree of regularity of 1(3^^, 1) from an explicit 
form of its Hilbert series. 

Recall that, if / is a homogeneous ideal of a polynomial ring R with ground field K, its Hilbert series 
is the series 

HSi{t)=Y,6:iuYvi{Rd/Id)t\ 

den 

where Rd denotes the K- vector space of homogeneous polynomials of degree d and Id denotes the K-vector 
space Rd n /. 

Proposition 1. The Hilbert series of the homogeneous ideal 1(3''^, 1) C ]K(o) [X] is 

where A{t) is the (p — 1) x {p — 1) matrix whose [i,j)-entry is (^^T*) (" fe "')^''- 
The proof of Proposition[T]is postponed to Section |33] 



Proof of TheoremU] By definition, the Hilbert series of a zero-dimensional homogeneous ideal is a poly- 
nomial of degree d^g —1- By Lemma|2l 1(5^'^, 1) has dimension 0. Thus, using Proposition[T] we deduce 
that: 



drcg(I(;?^l)) = l + deg 



The highest degree on each row of A{t) is reached on the diagonal. Thus deg(det A{t)) = ^'-^^ and a 
direct degree computation yields 

d.g(I(^^ D) = 1 + deg 

= D{p- 1) + (D - 2)n + 2. 

□ 



From Proposition [T] one can also deduce the degree of 1); this provides an alternate proof of 
Theorem 2.2]. 

Corollary 2. The degree of the ideal 1(3'^, 1) is 

DEG(I(ij'\ 1)) = (^^ 2 J j D^iD - . 

Proof By definition of the Hilbert series, the degree of the 0-dimensional homogeneous ideal 1(5^*^, 1) is 
equal to HSij^h By Proposition[T] direct computations show that HSijjh = det{A{l))DP {D ~ 

The determinant of the matrix A{1) can be evaluated by using Vandermonde's identity and a 
formula due to Harris-Tu (see e.g. ||22| Example 14.4.14, Example A. 9.4]). We deduce that det(A(l)) 



and hence HSi(s„,i)(l) = Q:^)^?^^ - 1)"-^- □ 
It remains to prove Proposition [T] This is done in the next subsections following several steps: 

• provide an explicit form of the Hilbert series of the ideal T>; this is actually already done in (TT\; we 
recall the statement of this result in Lemma|3] 

• deduce from it an explicit form of Hilbert series of the ideal 3(p„) using genericity properties satis- 
fied by the polynomials and properties of quasi-homogeneous ideals; this is done respectively in 
LemmalHand Section [X2l 

• deduce from it the Hilbert series associated to 1(3^'^, 1). 
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3.1 Auxiliary results 



We start by restating a special case of lfT2l Cor. 1]. 

Leimna 3 (|[T2 Corollary 1]). The Hilbert series of the ideal V C ]K[C/] is 

HSx)(i) = -j—r. ^— • 

Leiiuna 4. Foreach2 <l< np, does not divide in K{a)[U, X]/ 3i^f_iy 

Proof. According to [|29l Thm. 21 ll30l . the ring K(o) [U]/!) is a Cohen-Macaulay domain of KruU dimen- 
sion {n — 1 + p — {p — l)){p — 1) ~ n{p — 1) . Therefore, the ring K(a)[?7, X]/!) is also a Cohen-Macaulay 
domain, and has dimension np. 

Consider now the ideal (0i,...,0„p) C {K{a)[U]/V)[X]. According to LemmaEl the ideal 1(5^'^, 1) = 
(2? + {qi, . . . ,0„(p_i))) n K(a)[X] is zero-dimensional. Let -< denote a lexicographical monomial or- 
dering such that for all i,j,k, Uij >- Xk. Since the variables U can be expressed as functions of X 
- §§- ^ ^(p"))' we have LM^(2? + (oi, . . . , 0„p)) = {mj) + LM^(I(5'i>, 1)) which is zero- 
dimensional. Therefore, the ideal V + {q^, . . . ,g„p) C K(a)[[/, X] is zero-dimensional and hence so 
is {qi, . . . , g„p) C K(o)[C/, X]/!). Now suppose by contradiction that there exists £ such that divides 
in K(a)[?7, X]/ Let (q be the smallest integer satisfying this property. Since V is equidimensional 

and ye < 4,0^ does not divide in K(a)[C/, X]/ the ideal (g^, . . . ,0^„_i) C K{a)[U,X]/'D 
is equidimensional, has codimension Iq — 1, and thus has no embedded components by the unmixed- 
ness Theorem [[141 Corollary 18.14]. Since g^^ divides in the ring K{a)[U, X]/ {V + {gi, . . . , 
the ideal {g^, . . . C K{a)[U, X]/V has also codimension io — 1. Therefore the codimension of 

(fli5 • ■ • iQnp) ^ K{a)[U, X]/!) is strictly less than np, which leads to a contradiction since we have 
proved that the dimension of this ideal is 0. □ 



3.2 Quasi-homogeneous polynomials 

The degrees in the matrix whose entries are the variables Uij have to be balanced with D — 1, the degree 
of the partial derivatives. This is done by changing the gradation by putting a weight on the variables 
Uij, giving rise to quasi-homogeneous polynomials. This approach has been used in ||2T1 in the context 
of the Generalized MinRank Problem. A polynomial / G K[U, X] is said to be quasi-homogeneous if the 
following condition is satisfied (see e.g. [24, Definition 2. 11, page 120]): 

/(A'°""^wi,2, • ■ ■ , X^^^Up^n, Axi, . . . , Aa;*..) = A''/(ui,2, • . • ,Up^n,Xi, . . . ,Xk). 

The integer d is called the weight degree of / and denoted by wdeg(/). 

An ideal / C K[U, X] is called quasi-homogeneous if there exists a set of quasi-homogeneous gen- 
erators of /. We let K[U, X]'f^ denote the K- vector space of quasi-homogeneous polynomials of weight 
degree d, and /^"'^ denote the set K[U, n /. Ideals generated by quasi-homogeneous polynomials are 

positively graded, as shown in 11211 Proposition 1] that we restate below. 

Proposition 2 ( 11211 Proposition 1]). Let I C K[C/, X] be an ideal. Then the following statements are 
equivalent: 

• there exists a set of quasi-homogeneous generators of I; 

• the sets ij^^ are vector subspaces o/K[C/, X]''^\ and I = ^d"^- 

If / is a quasi-homogeneous ideal, then K[[/, X]/I is a graded algebra and hence its weighted Hilbert 
series wHS/(t) e ^[[i]] is well defined: 

wHS,(i) = ^dimK(K[C/,X]('"V4'"^)*'- 

d&i 

The following lemma and its proof are similar to 11211 Lemma 5] . 
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Lemma 5. The Hilbert series of 1(5^*^, 1) C K(o)[X] is equal to the weighted Hilbert series o/Jf^p^) ^ 
K{a)[X, U]. 

Proof. Let ^lex be a lex ordering on the variables of the polynomial ring K{a) [X, U] such that Xk ^lex Ui.j 
for all j. By inn Sec. 6.3, Prop. 9], HSi(5h4)(t) = HSLM^,^_^(i(5i',i))(t) and wHSa(p(„_i„ (i) = 
wHSLM-i|^ja(p(„_i)))(0- Since LM^,^Jmjj - fij) = utj and 3(p„) nK[X] = 1(3^'^, 1), we deduce that 

LM^,^p(j,„)) = /{u^.j-jU LM^,^p(p„)n]K(a)[X])) 
- <K,}ULM^,JI(5'\1))). 

Therefore, ) is isomorphic (as a graded K(a)-algebra) to . .. ^^^tS-^iI .i-, ■ 

Thus, HSlm^,^ji(5M))W = wHSLM^,^^(a(p„)) and hence HSi(5h^i)(<) = wHS^^^^^it). □ 

3.3 Proof of Proposition [1] 

We reuse Notations|2l 1(3^^, 1) = (2? + {qi, . . . , Qpn)) K{a)[X]. According to Lemma[3]and by putting 
a weight Z? — 1 on the variables U, the weighted Hilbert series of I? C K(a)[J7] is 

uc det A{t^-^) 



i(^-l)C'2')(l_iC-l)n(p-l)' 

Considering V as an ideal of K(o) [X, [/], we obtain 

wHSx,cK(a)[;7,X](i) = (1 - t)^ wHSpcK(a)[C/](0- 

If / C K(a)[/7, X] is a quasi-homogeneous ideal and if g is a quasi-homogeneous polynomial of 
weight degree d which does not divide in the quotient ring K(a)[C/, X]/I, then the Hilbert series of the 
ideal / + (g) is equal to (1 — f^) multiplied by the Hilbert series of / (see e.g. the proof of 11211 Thm 1] for 
more details). 

Notice that the polynomials fli, . . . are quasi-homogeneous of weight degree D — 1 (these 

polynomials have the form Uij — and the polynomials 0p(„_i)+i, • . • , 0p„ are quasi-homogeneous 
of weight degree D (these polynomials are j^i, • ■ • , fp). Since does not divide in ]K(ci)[f/, ^]/ 
(LemmalU), the Hilbert series of the ideal 3(p,i) C ]K(a)[X, U] is 

dctv4(i^-i)(l-<^)P(l-<^-i)p("-i) 

_ det A{t^-^) (1 ~ - 
~ ^(s-i)^^') (1-0" 
Finally, by Lemma|5] wHSaj^^j (i) — HSi(;jh 



4 The affine case 

The degree of regularity of a polynomial system is the highest degree reached during the computation of a 
Grobner basis with respect to the grevlex ordering with the algorithm. Therefore, it is a crucial indicator 
of the complexity of the Grobner basis computation. On the other hand, the complexity of the FGLM 
algorithm depends on the degree of the ideal I(F, 1) since this value is equal to dims (]K[X]/ I(F, 1)). 

In this section, we show that the bounds on the degree and the degree of regularity of the ideal 1) 
are also valid for (not necessarily homogeneous) polynomial families in K[X] under genericity assump- 
tions. 

Theorem 2. There exists a non-empty Zariski open subset 6" C K[X]^ such that, for any F in 6 H K[X]p, 

drcg(I(F,l)) < D{p-\) + {D -2)n + 2, 
DEG(I(F,1)) < 
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In the sequel, ]K[X]£i denotes {/ e K[X] \ dcg(/) = D}, and K[X]£)_hom denotes the homogeneous 
polynomials in K[X]£). In order to prove Theorem|2](the proof is postponed at the end of this section), we 
first need two technical lemmas. 

Lemma 6. There exists a non-empty Zariski open subset 6 C K[X]^ i^^^^ such that for all F'' S H 
IK[X]P LM^(I(F'\1)) = LM^(I(y'\l)). 

Proof. See e.g. lISTl Proof of Lemma 2] for a similar proof. □ 

Lemma 7. Let G = {gi, . . . , gm) be a polynomial family and let G'^ — {g'l, . . . , g^) denote the family of 
homogeneous components of highest degree ofG. If the dimension of the ideal (G'*) is 0, then DEG((G)) < 
DEG((G'')). 

Proof. Let -< be an admissible degree monomial ordering. Let LM^(ft,) denote the leading monomial of 
a polynomial h with respect to -<. Let m e LM^((G'')) be a monomial. Then there exist polynomials 
si,...,Sm such that LM^ (Sl^i ■^i9i) — Since ^ is a degree ordering, LM^ (X]"=i = 
Therefore LM^((G'')) C LM_^((G)). If the ideal (G'') is 0-dimensional, then so is (G) and hence 
DEG(LM^((G))) < DEG(LM^((G))). Since DEG(/) = DEG(LM^(/)), we obtain DEG((G)) < 
DEG((G'')). 

□ 

Proof of Theorem^ Let ^ be a degree monomial ordering, and F'' = . . . , ) € K[X]^ ^^^^ denote 
the homogeneous system where //* is the homogeneous component of highest degree of fi. By Lemma|6l 
there exists a non-empty Zariski subset C IK[X]c such that, forany F in fff^K[XY', LM^ (I(F'', 1)) 
LM^(I(5'^, 1)). By 1 13, Ch.9, §3, Prop.9], the Hilbert series (and thus the dimension, the degree, and the 
degree of regularity) of a homogeneous ideal is the same as that of its leading monomial ideal. Hence, by 
Lemma|2l 

dim(I(F'',l)) = dim(LM^(I(F'',l))) 
= dim(LM^(I(y'Sl))) 
= dim(I(5''\ 1)) = 0. 

Similarly, by Theorem[T] 

d„g(I(F^ 1)) = d,eg(I(;?^ i)) - d{p -\) + {d- 2)n + 2. 

The highest degree reached during the F5 Algorithm is upper bounded by the degree of regularity 
of the ideal generated by the homogeneous components of highest degree of the generators when this 
homogeneous ideal has dimension (see e.g. fS) and references therein). Therefore, the highest degree 
reached during the computation of a Grobner basis of I(F, 1) with the Algorithm with respect to a 
degree ordering is upper bounded by 

dreg < D{p -!) + (£>- 2)n + 2. 

The bound on the degree is obtained by Corollary |2] and Lemma |7] 

DEG(I(F,1)) < DEG(I(F'',1)) 

< DEG(LM^(I(?h^l))) 

□ 

5 Complexity 

In the sequel, w is a real number such that there exists an algorithm which computes the row echelon form 
of n X n matrix in 0(n") arithmetic operations (the best known value is a; w 2.376 by using Coppersmith- 
Winograd algorithm, see ||38l ). 
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Theorem 3. There exists a non-empty Zariski open subset C K[X]^, such that, for allY ^ 6 H K[X]^', 
the arithmetic complexity of computing a lexicographical Grobner basis o/I(F, 1) is upper bounded by 



D{p-1) + {D -2)n + 2 

Proof. According to Q H), the complexity of computing a Grobner basis with the Algorithm with 
respect to the grevlex ordering of a zero-dimensional ideal is upper bounded by 



O 



n + dreg 



dreg 

where dreg is the highest degree reached during the computation. In order to obtain a lexicographical 
Grobner basis, one can use the FGLM algorithm [,18] , Its complexity is O (nDEG(I(F, 1))'^) (better 
complexity bounds are known in specific cases, see ||T9l ). 

According to Theorem|2] there exists a non-empty Zariski open subset C K[X]^ such that, for all 

F in ^ n K[X]P, 

dreg(I(F,l)) < D{p-l) + {D-2)n + 2, 
DEG(I(F, 1)) < {'^zl)DP{D -ly-P. 

Therefore, for all F in ^ n ]K[X]p, the total complexity of computing a lexicographical Grobner basis of 
I(F,1): 

' 'D{p-l) + (L>- l)7^ + 2^ 



O 



D(p- 1) + (D - 2)n + 2 





(n- 




r- 












□ 



Corollary 3. If D — 2, then there exists a non-empty Zariski open subset G C K[X]2, such that for all 
F e ^ n the arithmetic complexity of computing a lexicographical Grobner basis o/I(F, 1) is 

upper bounded by 

o{{-^'^Wn2M--''" 



2p J \p- 1, 

Moreover, ifp is constant and D — 2, the arithmetic complexity is upper bounded by O (n'^P'^y 
Proof. This complexity is obtained by putting Z? = 2 in the formula from Theorem |3] □ 
In the sequel, the binary entropy function is denoted by h2'. 

Vx e [0,l],h2ix) = -a;log2(a;) - (1 - a;) log2(l - x). 

Corollary 4. Let D > 2andp € N be constant. There exists anon-empty Zariski open subset C K[X]^, 
such that, for allF € C\ K[X]p, the arithmetic complexity of computing a lexicographical Grobner basis 
o/I(F, 1) is upper bounded by 



Proof. Let a; be a real number in [0, 1]. Then by applying Stirling's Formula, we obtain that 

xnj \\/n 

Therefore, 
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Let C denote the constant D{p — 1) + 2. Then 

(D(p-l) + (D-l)n+2\ _ t(D-l)n+C\ _ ( ((D-l)n\\ 
\D(p-l) + (D-2)n+2) ^ \ n ) ^ ^ \\ n ) ) 

= O fj_2(^-l)''2(TrrT)"' 

The right summand in the complexity formula given in Theorem[3]is O (n^P{D — l)-^") when p and D are 
constants; this is upper bounded by 

O ( ^2(^-i)''^(^)" 



Let be the non-empty Zariski open subset defined in Theorem|3] For all F e n ]K[X]p, the arithmetic 
complexity of computing a grevlex Grobner basis of F is upper bounded by 

O (^^2^D-^Mo^h^^ = o [-jr^ii^ - l)e)"") 

^0({D~ i)(i+i/iog(£i-i))«c^) 

= o((i:'-i)3""), 

since D > 3 and lu < 2.376 with Coppersmith- Winograd algorithm. 

On the other hand the asymptotic complexity of the FGLM part of the solving process is 

O - 1)3") ^d(^{D- 1)3") , 

which is upper bounded by the complexity of the grevlex Grobner basis computation. □ 

The following corollary shows that the arithmetic complexity is polynomial in the number of critical 
points. 

Corollary 5. For D > 3, p > 2 and n > 2, There exists a non-empty Zariski open subset C K[X]^, 
such that, for F G D 'K[X]p, the arithmetic complexity of computing a lexicographical Grobner basis of 
I(F, 1) is upper bounded by 

O (dEG (I(F, 1))"MS^-'4)) ^ Q ^^^^ ^j^-p^ l))4 03-) . 

Proof. Let C K[X]^ be the non-empty Zariski open subset defined in Theorem|2] and F G ^ n K[X]^ 
be a polynomial family. First, notice that, since p>2 and n > 2, 

deg(i(f,i)) = (;:j)(i^-i)"-fi^p 

> n 

Therefore the complexity of the FGLM algorithm is upper bounded by 

O (n DEG (I(F, 1))^) < O (dEG (I(F, I))^) . 
The complexity of computing a grevlex Grobner basis of I(F, 1) is upper bounded by 

GREVLEX(p,n,i?) = O r(^(p-i)+(^^-i)«+2)") 

<o\nr). 

Notice that (^^") < {2DY'-^. By Stiriing's formula, there exists Co such that ^ < Cqg". Hence 
GREVLEX(p, n, D) = O ((2i5'e)"). 

Since Z) > 3 and ti < log(DEG(I(F, l)))/\og{D - 1), we obtain 



/ log(2er>) \ 

0((2£ie)"'^) < Om-E^"'^j 



, log(2eD) , 

< O DEG(I(F, l))i°s(i^-i)' 



The function D i-> is decreasing, and hence its maximum is reached for D = 3, and ^°^^^2) — 

4.03. °^ □ 
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Figure 1: Experiments in Magma measuring the arithmetic complexity (K = GF(65521)). 

Notice that in the complexity formula in Corollary |5] the exponent tends towards uj when 

D grows. Therefore, when D is large, the complexity of the grevlex Grobner basis computation is close to 
the cost of linear algebra O (DEG(I(F, 1))") . Also, we would like to point out that the bound in Corollary 

|5]is not sharp since the formula O (^("^j"^"*^) ^ for the complexity of the Fc, algorithm is pessimistic, and 
the majorations performed in the proof of Corollary 0are not tight. 

6 Experimental Results 

In this section, we report experimental results supporting the theoretical complexity results in the previous 
sections. Since our complexity results concern the arithmetic complexity, we run experiments where K is 
the finite field GF (65521) (Figure [U, so that the timings represent the arithmetic complexity. In that case, 
systems are chosen uniformly at random in GF(65521)[X]£i. 

We give experiments by using respectively the implementation of F4 and FGLM algorithms in the 
Magma Computer Algebra Software, and by using the F5 and FGLM implementations from the FGb 
package. 

All experiments were conducted on a 2.93 GHz Intel Xeon with 132 GB RAM. 

Interpretation of the results. Notice that the degree of regularity and the degree match exactly the 
bounds given in Theorem|2] In Figures [T]and|2] we can see a different behavior when D = 2 01 D = 3. In 
the case D = 2, since the complexity is polynomial in n (Corollary|3]l, the computations can be performed 
even when n is large (close to 20). Moreover, notice that for D = 2 or D = 3, there is a strong correlation 
between the degree of the ideal and the timings, showing that, in accordance with Corollary |5] this degree 
is a good indicator of the complexity. 

Also, in Figure |2] we give the proportion of non-zero entries in the multiplication matrices. This pro- 
portion plays an important role in the complexity of FGLM, since recent versions of FGLM take advantage 
of this sparsity 1 19 |. We can notice that the sparsity of the multiplication matrices increases as D grows. 

Numerical estimates of the complexity. Corollary|5]states that the complexity of the grevlex Grobner 
basis computation is upper bounded by O (DEG(I(F, 1))* °'^") when D > 3, p > 2, n > 2. However, the 
value 4.03 is not sharp. In Figure|3] we report numerical values of the ratio log ("^,'J'°*') / log (DEG(I(F, 1))) 
which show the difference between 4.03 and experimental values. 

Notice that all ratios are smaller than 4.03, as predicted by Corollary |5] Experimentally, the ratio 
decreases and tends towards 1 when D grows, in accordance with the complexity formula 

(log(2eI?) \ 
DEG (I(F,l))'°«("-i)"J 
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Figure 2: Timings using the FGb library and K = GF(65521). 
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for the grevlex Grobner basis computation. Also, when D > 3, the worst ratio seems to be reached when 
p = n— 1, D = 3 and n grows, and experiments in Figure|3]tend to show that it is bounded from above by 
2.28. 
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